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ONTARIO MINISTRY OF THE ENVIRONMENT 
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ERRATA 



"Sound Attenuation in Lined Ducts, Simplified Design Pro- 
cedures", August 198 1, J. Buur, C, A, Krajewski. 

1. Line 14 from the bottom on page 1 (just under Figure 
1.1) should read i 
"The attenuation D in dB .»•..*•..,.•* " 

2® Equations (§} , (7), (8) and (9) on page 6 are expressed 

in G.G.S. units . M.K.S, form of the equations is given 
on page 257 Of the Reference 15. 

1, Lines 3, 4 and 5 on page 9 should reads 

"The resulting constriction of the flow causes an 
increase in acoustic impedance." 

4. The paragraph under Figure 2.3. on page 9 should read: 

6 

"In such case, according to ingard , the resistivity is 

increased by R f 6 , where the "end correction" 

6- R_ s f (p) is found according to the open area ratio p 

of the perforation, using Figure 2.4." 

The next sentence, as well as the formula shown beside 

Fig, 2.4 should be deleted. 






c, * 



2 - 



In aquation (13) on page 10 the term r should be replaced 
by R f . 



6. Term D. in equation (23) and in line 15 on page 13 should 
be replaced by D. 

7 . In Figures 3.3, 3.4. and 3.5. the input parameter 

5 
R-l, 48 x 10 mks rayls/m should be replaced by equation 

(33) from page 17. 

8« Line 4 on page 23 should read: 

"c sound velocity (m/s) " 

9, Line 6 on page 23 should read: 

"D, attenuation per duct length (dB) " 
equal the duct height. 

10. Line 9 on page 23 should read: 

"k wave number = u (radians/m) " 
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SUMMARY 

This is an attempt to review methods for the determination 
of acoustic impedance of various absorbent lining designs 
and the available procedures for calculating duct attenua- 
tion, A comparison of the available methods and the suggested 
design procedures are included in the body of the report. 

The draft report is the work of Jacob Burr, a foreign student 

employed by the Ontario Ministry of the Environment, Toronto, 
under the International Association for the Exchange of 
Students for Technical Experience (IASTE) exchange program. 
At the time of his short 1981 (six weeks) work term with the 
Noise Pollution Control Section he had completed three years 
Q'f training in electronic engineering at the Technical 
University of Denmark in Copenhagen. 

The project was assigned, and work progress was supervised 
by C. A. Krajewski, P.Eng. , Senior Acoustical Engineering 
Specialist in the Noise Pollution Control Section, Pollution 
Control Branch, 
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1. INTRODUCTION 

The problem of noise propagated through a duct which 
allows the passage of air or other fluids is usually associated 
with air handling or combustion equipment, such as? boilers, 
fans, combustion engines, etc. An effective method of 
reducing such noise is to line one or more of the inner 
walls of the duct with sound absorbing material, see Figure 1.1, 




Figure 1.1. 



V" 



The attenuation D in db afforded by such duct is a 
rather complex function of the duct geometry (i.e. the free 
height h and length 1), the frequency of the sound, the 
temperature of the flowing medium, and the acoustic impedance 
of the absorptive lining material. 

A number of approximations to the exact theory for 
calculating this attenuation exists, since the attenuation 
equation can only be solved analytically for simple cases. 



An analysis of acoustic effectiveness for the lined 
duct involves two uncorrelated steps: 

«- Determination of the acoustic impedance Z for the 
absorbing wall(s). 

- Calculation of duct attenuation based on its 

geometry and known value of the acoustic impedance 
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As no comprehensive summary of design procedures for 
various types of ducts can be found in the available literature, 
an attempt was made in this paper to: 

(i) review theory related to determination of the impedance 
for various absorbent lining designs; 

(ii) review basic theory of sound attenuation in ducts; 

(iii) compare various approximate calculation methods; 

(iv) provide machine format for a practical design assess- 
ment procedure. 

This project is limited to theoretical analysis as no 
duct equipment or measurement facilities are available for 
testing. It is suggested that compilation of experimental 
results from various sources on attenuation for lined ducts 
be an extension of the present study. 



2. DETERMINATION OF IMPEDANCE FOR THE DUCT LINING 



2,1. Acoustic Impedance 



of air, is: 



Z m p «c m 1.20kg/m 3 x 343 m/s = 413 kg/m 2 s 
00 o 

(at t = 20 C) 



w 






The specific acoustic impedance Z of an absorbing wall 
is defined as the ratio of sound pressure p to the particle • 

velocity u at the same point of the medium, 

Z - -£- (kg/ 2) or (mks ray Is) (1) 

u si J 



Z = R + jX 

The impedance is a complex variable and consists of the 
resistance R - Re Z and the reactance X = Im Z . The acoustic 
impedance of air Z, often called the characteristic impedance 
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It is convenient at times to express the impedance as 



the normalized impedance 2 



n 



n 



o 



p c 





(2) 



2,2, Lining Design 

The design of the lining is an important factor in 
achieving desired duct attenuation. Type of lining is 
selected to provide the optimum impedance at the frequency 
where greatest attenuation is required, However, other 
parameters such as temperature, pressure drop and dust load 
of the passing media must also be considered in the selection 
process. Figure 2.1. shows various commonly used lining 
designs. 
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Figure 2.1 



For a homogeneous bulk of absorptive material in front 
of a rigid wall (see configuration jj , the resistance is 
determined solely by the material, and the reactance is nega- 
tive in all practical cases, as the air in the layer acts 
very much like a spring because of the rigid wall. 
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In the case of a thin layer of absorptive material at a ^ 
distance d from the wall (see configuration ii) /the air 
cavity adds to the imaginary part of the impedance. 
To avoid propagation of sound behind the absorbing layer, 
the airspace is often subdivided into cells by partitions. 

By applying a thin sheet of perforated (hard) facing to 
absorptive material (see configuration iii) ,both the resistivity 
and the resistance changes. Most important, the mass of air 
in the holes of the perforated facing results in a positive 
contribution to the reactance, thus opening the possibility 

of tuning the system to a required resonance frequency. 

-• 

2.3 Parameters Describing Acoustic Behaviour of Absorptive 

Materials . 
Porous material used for the purpose of attenuating sound is 
normally composed of fibers (mineral wool, glass) or foam. 
To describe the acoustic behaviour of absorptive materials 
the following parameters are used: 

a) Flow Resistivity R f defined as the ratio of sound 

pressure change through the material 

to the particle velocity u times 

thickness 1. ^ 

AP 3 

R.e s (kg/sm ) or (mks rayls/m) * 

u»l 

b) Porosity P defined as the ratio of the volume 



of voids V to the total volume V 
a m 



P = 



V 
m 



In the case of fibrous material the porosity is 
expressed as: 

P = 1 - 



p m 



P 



where p is the density of the 
m 

material and p. is the density of 
the fibres. 



w 
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c) Structure Factor m accounts for the increase in 

inertia of the gas in the material 
due to the inner structure of the 

bulk. There seems to be a way of 

1 5 
measuring m, but Beranek * has 

given an approximate relation 

between structure factor and porosity 

of a homogenous material. See 

Figure 2.2. 

Structure factor m 
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Figure 2.2. 

2.4 Impedance of Absorbing Lining, 

The acoustic impedance of an absorbing lining may be 

measured or approached theoretically. 

9 
Kurze gives a method for measuring the complex impedance, 

12 
and Dawern ' presents an impressive collection 

of data measured on a duct silencer system consisting 

of a blanket of porous material with perforated facing, 

backed by an airfilled cavity, (Figure 2. I.e.) 

A theoretical analysis of acoustic impedance for dissipative 



lining materials has been conducted by Bolt, 1947 , Ingard 

8 



3 4 

and Bolt 1951 , Callaway and Ramer 1952 , Ingard 



6 7 

1954 , Kurze 1965 , and Velizhanina 1965". It generally 

involves the solution of the wave equation for a layer 

of absorbing material of thickness d backed by a rigid 

wall: 



Z = 



(3) 



tanhyd 
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where Z is the characteristic impedance of the material and y is 



the propagation constant of the material, 
These coefficients are given explicitly b; 
resistance r , porosity P and structure factor m as defined earlier: 



8 
These coefficients are given explicitly by Velizhanina , using flow 




( 4 ) 



and 



T = J "l/lT 2 I "f- P - 3 ~ ) (5) 



where q is the specific heat ratio at constant pressure and 

constant volume, and K is the compressibility of air. 

For isothermal conditions in the material q=1.4 and for 

adiabatic condition q-1 . 

At low frequencies the conditions are likely to be isothermal. 

Through a number of measurements on a wide range of absorptive 

1 
materials Delany and Bazley have empirically determined 

power relations for the characteristic impedance and the pro- 
pagation coefficient which are both functions of the flow 
resistivity R f : 



Z Q = R+jX , y = a+ j, 



f -0.75 
R = p c + 9 .08 p c C-S — ) (6) 

R f 



„ -0.73 
X = - 11.9 pc (-|— ) (7) 



-0.59 
I 
1 



a - 10.3-^- (-§-) 






* 
*> 



Thus, there are two different approaches to the 

determination of value Z for an absorptive lining, 

using either the analytically or the empirically derived 

functions for Z and Y • 
o 



2.5 Thick Layer of Absorbent J&terial, 

In the case of a thick layer of porous material backed 

by a rigid wall, as shown in Figure 2.1. a. , equation (3) 

may be applied directly. 



2.6 Thin Layer of Absorbent toterial Backed by Air-cavity. 

This type of wall design, shown in Figure 2.1.b, has, 

8 

accordingly to Velizhanina , the acoustic impedance: 



Z - 



TTq 



{ Pql * d) + j 



m 



pm 



pc 



3 P q 



T-^ 1 + d > " a^ql+d) 



(10) 



The formula is the result of expanding the tanh in a 
series and retaining only the two first terms. The 
last part of the reactance represents the stiffness of 
the air in the cavity between material and wall, similar 
to the formula for a Helmholtz resonator , (Pql + d) 
being the depth of the resonator. The first part of 
the reactance is the mass of air ( P*l) moving in 
the layer of material. The "torturous path" the air has 
to make through the layer is taken into account by 
multiplying the first component of the imaginary part by 
the structure factor m. At low frequencies this component 
will normally be small compared to the stiffness of 
air. 



Ingard and Bolt give a simpler approximation of the 
impedance of the same wall configuration ,as follows: 
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Z - rl + j 



.p ml - 



P c 



tan(-^-d) 



(11) 



This equation is derived on the assumption that the 
porous layer is very thin compared to the wavelength of 
sound, i.e. kl<<l .in such case the resistivity of the 
material is simply rl. It is recognized that the 
stiffness of air in the cavity is here expressed by 

the equation used for the one-dimensional case of a 
piston radiating into a narrow tube 



P c 



tan ( oj/c) d 



This holds only if the partitions dividing the cavity into 
cells are close enough together compared to the wavelength, 
so that wave motion may be assumed in one direction only. 
This is likely to be the case when spacing b< /4 
If b is larger or no partitions are present, the term becomes 



P C 



tan( -7r-d cose )cose 



(12) 



where 6 is the angle of incidence of plane wave 
(normal incidence: 8=0^ and cos8=1 ). 

This expression may only be used at low frequencies as 

shown by comparing the expressions below ( See Appendix A.I.) 



and 



to ' 



ta 
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2.7 Perforated Facing 

Perforated facing applied to an absorbing wall distorts 
the flow in a small region around the perforation. The 
resulting constriction of the flow causes an increase 

in acoustic mass reactance. 

If the facing is in close contact with the porous 
material i.e. the airspace between facing and material 

§ is smaller than the diameter 2a of the holes 
(see Figure 2.3.i±)and iii) the effect of the facing is 
not only mass reactance but a resistance as well. 



i\ & 




ii) V/////////////////f/A 




iii ) 




s 2a 



Figure 2,3, 



In such case, according to Ingard r the resistivity is 

increased by ' \ where the "end correction* 1 ' | is 

P 
found from Figure 2,4 • according to the percentage open 

area of the perforation p , The resistivity may be 

even larger than the resistivity due to the material. 



Up) 



1 I 1 ,i 1 1 ', 1 


• r I i , | 
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1 . i i ! i | i i 


1 — r~*'"^-T-— j 1 1 



6 - S»f ( P ) 

where s is the dis- 
tance between holes 



Figure 2,4. 
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The mass of air in the holes results in the reactance 
(t + (1+m) 6)/p to be added to the impedance of 
the wall without facing Z fe , as expressed by equations (3), 
(10) and (11) : 



Z = Z b + rs/p* + j wp (t + (1+m) 6) / p| (13) 



If the facing is not in close contact with the absorbent 
layer, the resistance contribution to the impedance is 
negligible, but reactance remains the same. 

8 
Velizhanina uses a slightly different "end correction" 6 

v 

given by: 

f - 0.85 a ( 1 - 2.50 a/s) (MJ 

v 

Thus, the impedance is expressed by: 

Z = Z y + r 6 y /p* + jup ( m + 1 ) « v / P° (15) 
where Z is found in equation (10) 

A viscous resistance from the holes is also present, 

although it is negible in most practical cases. According 

2 
to Bolt the value of viscous resistance is approximately: 

R y ■ 4 t/aJr P q b f (16) 

where q. is the coefficient of viscosity of air. 
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2.8 Other Types of Lining, 

7 
Kurze qives feres other examples of absorbent walls 



///// 



////- / 



// 



SlJi 



Figure 2,5, 

(i) The thin layer of absorbent material is located halfway 
between the wall and the end of partitions as shown in 
Figure 2.5. For this type of lining the impedance is 
given by: 



1 - cot kd/2 - jrl cot kd/2 
rl - j 2 cot kd/2 



(17) 



(ii) A set of Helmhoitz resonators with resonant frequency 



res 



, as shown in Figure 2 . 6 . 




For this type of lining the impedance is given by. 



1 .R-b 2 + j p C 



res 



res 



res. 



(18) 



where 



R - — -— p -Jz wv 

it a 



JL + 2 ( i - Utf. , 



n 



3. 
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CALCULATION OF DUCT ATTENUATION 



The theory of the transmission of sound inside ducts 

1 5 

has been studied by Morse, 1939 , Cremer, 1953' and 

7 § 
Kurze, 1965 and 1969 . All three considered the case 

of a duct of rectangular cross-section with one ab- 
sorbing wall and the other three walls perfectly rigid 
The derived formulas also apply to the case of two 
equally lined walls facing each other on basis of the 
mirror image principle, provided the dimension h, the 
free duct height, is defined as in Figure 3.1, 



|BMffiOffi& 




1 - 


is 


mmmi 





2h 



Figure 3, 1 

The attenuation of treated duct is usually expressed 
as in dB per unit length equal to free duct height h 

3.1 Morse and Cremer Method 

Morse solves the wave equation for two dimensional 

case (x, y coordinates) and the simplest type of wave: 

p h P Q cosh I (iry/h ) (k-jp ) exp j u t - (u/c) (a+jx )x (19) 



defining the transmission parameter (o+jr ) and 

the distribution parameter (u+jic) as follows: 



(a+jx) 2 ■ f{p/n) + j(«/n) ] 2 - 1 

= l/ (p + j'K 



(20) 



n (cj+JT ) 



,2 2 

J ~ Hi 
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where n is the frequency parameter; 

2hf 



(21) 



Applying a set of boundary conditions to eq. (19) , 
Morse determines the relationship between the distribution 
parameter ( y+js) and the acoustic impedance Z of the 
absorbing wall: 



pen 



cot h [- jTi(y + j(c)] ( 22 j 

y + j K 



The attenuation in dB/m is determined by: 



D h * 8.686irri<J = 8.686* Re ( J(u + j*) 2 - n 2 ) ( 23 J 



Since it is not possible to find (y+j»c) analytically 
from equation (22) , Morse introduced nomograms which 
show the distribution parameters for a given impedance 
| z|, phase angle <f> and frequency parameter n . 
The nomograms, however, are not suitable for design 
purposes, since D, still has to be found by applying 
equation (20) .Therefore ,Cremer presented graphs which 
yield the value of a for a set of values for n , 

p-^rT and #« 

Furthermore, Cremer shows that types of waves of higher 
order than the one considered by Morse will be attenu- 
ated more effectively than calculated. On the other 
hand, high frequency waves will tend to form "rays" 
close to the center of the tube. This effect causes the 
attenuation D to decrease proportional to n and even 



_ 2 

fj . This effect should be considered if 



-H- 



P c 



£> (n + l/2)<< 0.3 



(24) 



where n is the mode of the sound wave , 



Further , Cremer states that while 



p c 1 
— n * T 



(Z5) 



the problem may be solved as a one-dimensional case 
and thus equation (22) is simplified to: 



lv * 4*3 = I — 



•pc 



(26) 



3 . 2 Beranek Method 

For the study of duct silencers, using various 

7 
absorptive lining configurations, Kurze has slightly 

5 

altered Cremer 1 s graphs to show actual D, values in 

dB's instead of the parameter B • These graphs 

15 
have been used by Beranek to develop a procedure 

for calculating the attenuation of a lined duct, The 

procedure is outlined in the following steps: 



1. calculation of the impedance Z at n 



2, 

3. 



calculation of modulus |z| and argument 9 

determination of the appropriate value of D, 

III n 

1 ' ' and 6 to the graph (see Figure 3.2) 



ing 



pen 




Figure 3.2 - 

One of a number of graphs drawn by Kurze to obtain D, 

as a function of the lining impedance. 



- 15 - 

This method represents a fair approximation to exact theory 
but it is rather cumbersome for practical applications. 
By computerizing the first two steps, D, can be evaluated 
for a number of data sets in order to compare the methods 
discussed in the following two sections, 



3.3 Kurze Method 

Kurze derives a formula based on an assumption that the 
lining is of periodic structure. Considering the sound 
waves only at a number of discrete points and applying 
complicated matrix calculation the following equation is 
derived: 



v = 



D ( r + s 
v v 



4nZ 



2nZ + j kh pc 



) - 



(27) 



u n^ / k-h \ 2 t y h \ 2 
where D = ( — — ) - 2 + (*^^) 



and r , s are the coefficients of the Tschebyschef f 
function associated with this equation, and n is the 
number of discrete points along the dimension h. Letting 

n * » , the equation (27) approaches the equation derived 

15 
by Morse and Cremer . 



According to Kurze, the equation (27) , for n = 2, yields 
results within 5% of those obtained using Beranek's 
graphical method. The complex propagation constant given 
by Kurze is shown in the following equation: 



Y = 




( 2 ) 2 



1 + 



1 + 



47 
jkhp c 



(1+ jkhpc' 



(28) 
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+ 
The _ sign in the equation (28) is chosen so as to yield 

the smallest attenuation D, , which is given by; 



D h = 8.686 Ref y ft J [dB] 



1 
Substituting H = jkhpc— s~ * duct attenuation is expressed 

by: 



0. s 8,686 |il| E - (kh) 2 I (29) 

where E m 16 + 8H | V256 + 128H + 32H 2 j 3Q j 

4 + H 



Equation (30) is in a suitable forza for machine analysis. 



3 . 4 Mechel Method 

A formula similar to equation (29) but containing different 
numerical coefficients is given by Mechel as follows: 



^E - (kh) 2 



D. - 8 . 686 Re { m t - | kh j j (31) 

where 



E . 105 + 45H I yil025 + 5250H + 1650H 2 (32 ) 

20 + 2H 



Dividing E by 6,5, it becomes clear that Mechel's 
coefficients are close to those shown in equation (30) 
except for denominator 20 + 2H; 



E s 16.2 + 6.9H j 7260.9 + 124, 3H + 39. 1H 2 
3.1 * 0.3H 
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This indicates that the formula is derived in a way similar 
to Kurze's. Unfortunately, details O'f Mechel's theoreti- 
cal derivations were not published and it is, therefore, 

difficult to determine limitations of his method. Accord- 
ing to Mechel, the calculated values are always within 
5% of the exact solution. 



4 , CO MPARISON J3F_ METHODS 

For the purpose of comparison among the above discussed 
methods, the following simplified impedance equation 

was used: 



Z = 3 pc - j 46 pc/fd (33) 



i.e. the same as given by Beranek 



Attenuation values, calculated by all three methods for 
the ratio d/h of 0.25, 1.0 and 4,0 are shown in 

Figures 3.3, 3.4, and 3.5. The results obtained using 
Kurze's formula and simplified Mechel * s formula are 
almost identical in the low and medium frequency range. 
Slight discrepancies appear above the frequency of 
1 Khz, with Kurze* s formula predicting approximately 
0.08 dB (per unit length equal to the height of the duct) 
lower values than Mechel 's simplified formula. 
Attenuation values obtained using Beranek 1 s nomogram 
method are, in the low frequency region up to 500 Hz, 
generally lower than those predicted using the other two 
methods by 0.15 to 0.2 dB. This seems to be compensated 
by higher attenuation values in the frequency region 
above 500 Hz. The closest agreement among all three 
methods is noted for higher ratio of d/h, and above 
the frequency of 500 Em, 




00 
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I— » 

to 
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3. CYLINDRICAL DUCTS 

Morse approaches the problem of ducts with circular 
cross-section and gives the solution in the form: 

Z _ _ A - jB 



pen -j (C-jD) (y+j K ) 



134) 



where (p+j k ) is defined by equation (26) and A, B, C, 
and D are all series of Bessel functions . 

11 
Mechel derives a formula quite similar to his own 

formulas for the rectangular duct. Only the quantity h 

was changed from duct height to radius, and the numerical 

coefficients are different. It is presented in the 

following form: 



D h » 8.686 Re I Je- <kh) 2 } (35) 



where 



96 + 36H ±19216 + 2304H + 912H 2 

E - — ■ 

12 + H 



A comparison between results of calculation {using the 
equation 35) and experimental results of measurements 
presented by Mechel, show rather good agreement. 



6, OPTIMAL WALL IMPEDANCE 

1 
According to Cremer , the impedance 

Z = { 0.91 - j 0.76 ) pen (36) 

yields the highest attenuation. 

Since this attenuation is rather difficult to achieve 
over the entire frequency range, the design of a duct 
silencer normally involves selection of a lining which 
will provide the optimum impedance at a given resonance 
frequency . 
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The highest attenuation values will be achieved with the 

5 
lowest acceptable value of duct height h. Cremer also 

states that attenuation is greater the "closer the 
absorbent lining is to the sound propagation" . A rectan- 
gular duct is, therefore preferred to a round one. 



7. SUGGESTED DESIGN PROCEDURE 

7.1. Calculation of acoustic impedance for the selected - 

duct design: 

case of: * 

(i) Homogeneous bulk of absorptive material in 
front of a rigid wall. (Figure 2.1. a), 
Use equations (3), (6), (7), (8), and (9). 

(ii) Thin layer of absorptive material at a dis- 
tance d from a rigid wall (Figure 2,l,b.), 

Use equation (11) 

(iii) Thin sheet of perforated facing applied to 
absorptive material (Figure 2. I.e.), 



» 



Use equation (13) 

(iv) Absorbent wall (Figure 2.5.), 

Use equation (17) f 

(v) Set of Helmholtz resonators (Figure 2.6.), 
Use equation (18) 

7.2 Calculation of duct attenuation: 

(i) Rectangular ducts, 

Use equations (28), (29), and (30). 
(ii) Cylindrical ducts, 

Use equations (35) and (36) 






* 
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LIST OF VARIABLES UNIT 

a radius of holes in perforated facing 
b space between partitions 
c sound velocity 
d thickness of absorbing wall 



attenuation per duct length = D— s- 



D h 

f frequency 

h height of duct 

k wave number = - 



( 


m ) 


( 


m ) 


( 


m ) 


( 


m ) 


( 


dB/m ) 


( 


s ) 


( 


m ) 


( 


m ) 



c 

K compressibility of air 

1 thickness of absorbing blanket ( m ) 

m structure factor 

2 

p sound pressure ( N/m ) 

R- flow resistivity ( mks rayls/m ) 

r flow resistance { mks rayls ) 

t r temperature ( °C ) 

u particle velocity ( m/s ) 

I X acoustic reactance ( mks rayls } 

Z 



z normalized impedance = ■ 



P-c 



2 
Z specific acoustic impedance ( kg/m s ) 

2 

Z characteristic impedance ( kg/m s ) 

o 

n frequency parameter = — — ■ 
Y propagation constant 

m+Jk distribution parameter 

3 
P density of air ( kg/m ) 

<?+J T transmission parameter 

A wavelength ( m ) 

x co-ordinate parallel to duct axis 
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